Resonant and rolling droplet 
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When an oil droplet is placed on a quiescent oil bath, it eventually collapses into the bath due 
to gravity. The resulting coalescence may be eliminated when the bath is vertically vibrated. The 
droplet bounces periodically on the bath, and the air layer between the droplet and the bath is re- 
plenished at each bounce. This sustained bouncing motion is achieved when the forcing acceleration 
is higher than a threshold value. When the droplet has a sufficiently low viscosity, it significantly 
deforms : spherical harmonic Y/^ modes are excited, resulting in resonant effects on the threshold 
acceleration curve. Indeed, a lower acceleration is needed when £ modes with m = are excited. 
Modes m 7^ are found to decrease the bouncing ability of the droplet. When the mode £ = 2 
and m = 1 is excited, the droplet rolls on the vibrated surface without touching it, leading to a 
new self-propulsion mode. 

PACS numbers: 45.55.Dz,82.70.Uv,83.80.Qr 



Bouncing droplets have been first studied by Couder 
a/ With a radius R < 0.5 mm, they may gener- 
ate local waves on a 50 cSt bath and use them to move 
horizontally and interact with their surroundings : they 
detect submarine obstacles, they orbit together or form 
lattices, they are diffracted when they pass through a 
slit [21, 0, 0, S @1 • Bouncing droplets exhibit a wide vari- 
ety of motions and interactions that present some strong 
analogies with quantum physics, astronomy and statisti- 
cal physics. 

In contrast with these previous studied, we choose a 
bath viscosity of 1000 cSt, a droplet viscosity between 
1.5 and 100 cSt and a droplet radius R between 0.76 
and 0.93 mm. With a bath viscosity at least 10 times 
larger as the droplet viscosity, the bath deformations are 
inhibited (capillary waves are fully damped) while the 
harmonic deformations of the droplet are enhanced. As 
shown recently in [8], the droplet deformation ensures 
its bouncing ability. Various modes of deformation may 
be excited as depicted in FiglH Each picture has been 
constructed from an experimental snapshot of the droplet 
(the left side of each picture) and the calculated 3D spher- 
ical harmonic (the right side). Those modes are analo- 
gous to the natural modes of deformation introduced by 
Rayleigh [10] and may be expressed in terms of spherical 
harmonics Yf^. The simplest mode that may be used for 
bouncing is the mode Y2 [8] . We will show that droplets 
can use non-axisymetric mode Y2 to move horizontally, 
or more precisely to roll over the bath. This new mode 
of self-propulsion drastically contrasts with the bouncing 
walker mode described in Q since the bath deformations 
are not necessary for the propulsion. 

On a bath vertically vibrated according to a sinusoidal 
motion Asin(27r/t), the periodic bouncing only occurs 
when the reduced maximal vertical acceleration of the 
bath T = Air'^Af'^/g is higher than a threshold value Tth^ 
where g is the acceleration of gravity. Formally, Tth de- 
pends on the forcing frequency / [8], the droplet radius 
R [7], and the physical parameters of the liquids (den- 
sity p, viscosity surface tension a). We measure the 
threshold Tth as a function of the forcing frequency for 



FIG. 1: Various deformation modes of a bouncing 
droplet (ly = 1.5 cSt, R = 0.765 mm) observed with 
a high-speed camera for various forcing frequencies. 

The first line (resp. 2nd, 3rd) displays a mode Y2 (resp. 1^3°, 
Y4) with m = (axisymmetry). The forcing frequency is 
50 Hz (resp. 160 Hz, 275 Hz) and the reduced acceleration 
r = 0.3 (resp. 2, 6). First and second columns represent snap- 
shots at two different times of the oscillation. The spherical 
harmonic solution (on the right of each picture) is superposed 
to the experimental pictures (on the left of each picture). 



various droplet sizes and viscosities. For each frequency, 
we place a droplet on the bath in a bouncing configu- 
ration, i.e. F > Tth' The forcing acceleration is then 
decreased progressively. When the threshold is reached, 
the droplet cannot sustain anymore the periodic bounc- 
ing and quickly coalesces with the bath. Since the droplet 
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and the bath are made from different hquids, a coales- 
cence event locally contaminates the bath. Therefore, 
droplets always need to be placed at different locations 
on the bath. 



102 



101 



10° 



10- 



■ • 1.5 cSt 
- • 10 est 

■ • 100 est • , 


.V* 

^ V : 






♦ 

yO 
^3 














Roller 







50 100 150 200 250 300 350 
f(Hz) 

FIG. 2: Evolution of the bouncing threshold acceler- 
ation Tth with respect to the forcing frequency. The 

red, green and blue bullets correspond to droplets with viscos- 
ity ly =1.5, 10 and 100 cSt respectively, and with a constant 
radius R = 0.765 mm. Depending on the forcing frequency, 
various modes are observed, that may be related to spherical 
harmonics . 

In FiglSl the threshold acceleration Tth is represented 
as a function of /, for droplets with viscosities u =1.5, 
10 and 100 cSt. The threshold acceleration for the 100 
cSt droplet increases with the frequency in a monotone 
way. By opposition, the 1.5 cSt curve is characterized by 
regularly spaced local minima that, in some way, corre- 
spond to a resonance of the system : a minimal energy 
supply is required to sustain the periodic bouncing mo- 
tion. The first minimum, at / = 50 Hz, corresponds to 
Tth =0.25, a value significantly less than the 1 g mini- 
mal threshold required by inelastic bouncing objects on 
a vibrated plate. Modes m = and i =2, 3 and 4 are ob- 
served for forcing frequencies corresponding to minima in 
threshold acceleration curve. The higher the frequency, 
the higher-order the excited mode, and the higher the 
corresponding threshold acceleration. We note that the 
asymmetric mode Y2 occurs at a local maximum of the 
threshold curve. 

The bouncing droplet may be considered as an oscil- 
lating system analogous to the damped driven harmonic 
oscillator : surface tension is the restoring force and vis- 
cosity the damping process. The dimensionless ratio be- 
tween both is the Ohnesorge number Oh = v^j \J cfR^ 
which is equal to 0.012, 0.078 and 0.775 for a droplet vis- 
cosity of 1.5, 10 and 100 cSt respectively. When O/i <C 1, 
the viscous damping is negligible and resonance is im- 



portant. Damping increases as Oh gets closer to 1. As 
seen in FiglSl the bouncing droplet reacts to an increased 
damping in the same way as the harmonic oscillator : {%) 
the resonance frequency slightly decreases, shifting the 
whole threshold curve to the left, {%%) the required in- 
put Vth increases, and (m) extrema tend to disappear. 
At 100 cSt {Oh = 0.775), the damping is fully active 
and the threshold curve increases monotonically with the 
frequency : no more resonance is observed. This high- 
viscosity behaviour has already been observed by Couder 
et at [l|] for 500 cSt droplets bouncing on a 500 cSt bath. 
Those authors proposed to model the threshold curve 
by Tth = 1 + OLp where a depends on the droplet size 
among others. This equation is obtained by only consid- 
ering the motion of the mass center of the droplet and 
the squeezing of the air film between the droplet and the 
bath, without considering the droplet deformation. This 
model only fits high viscous regimes in contrast with the 
model developed in [8]. 
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FIG. 3: Dimensionless amplitude Athfc/o a func- 
tion of the reduced forcing frequency f j fc-> where Ath 
is the amplitude of the bath vertical motion corresponding to 
the reduced acceleration Vth and /c = \J M/af the capillary 
frequency of a droplet of mass M = 47T/3pR^. Blue, green 
and red bullets correspond to droplet radius R = 0.765 mm, 
0.812 mm and 0.931 mm respectively. The droplet viscos- 
ity is = 1.5 cSt. Modes ^ = 2, ^ = 3 and ^ = 4 are 
observed in the black, yellow and red ranges of frequency. 
Boundaries between those zones correspond to maxima in the 
threshold curve. Moreover, those particular frequencies may 
be obtained by multiplying the Rayleigh natural frequencies 
(EqU]) by a factor 1.15. Arrows indicate the Rayleigh fre- 
quencies, which cannot be directly related to inflections in 
the threshold curve. 

In 1879, Lord Rayleigh [T^ described natural oscilla- 
tions of an inviscid droplet. Since those oscillations are 
due to surface tension, na tural frequencies scale as the 
capillary frequency fc = ^JajM^ where M = Air/SpR^ is 
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the droplet mass. More exactly, the dispersion relation 
prescribes the natural "Rayleigh" frequency fn related 
to a £-mode: 




The function F{£) may vary (frequencies are shifted by a 
multiplicative factor), depending on the way the droplet 
is excited [11, 12]. For free oscillations, Eq.([Tj) is degen- 
erated according to the m parameter. In Fig. [3l thresh- 
old data obtained for various droplet sizes collapse on 
a single curve by using the Rayleigh scaling. Moreover, 
at the bouncing threshold, the vertical force resulting 
from the droplet deformation exactly balances the grav- 
ity. One may define a characteristic length L = g/f^ 
corresponding to the free fall distance during the capil- 
lary time 1/ fc As shown in Fig. [3l the threshold ampli- 
tude Ath = Tthg/i^Tr'^f'^) scales as the length L, what- 
ever the droplet size. The minimum value of Athfc/d 
does not vary significantly with the mode index £. In 
Fig. [31 the natural frequencies defined in Eq.([T]) and rep- 
resented by arrows do correspond neither to the minima 
in threshold, nor to the maxima. However, these fre- 
quencies multiplied by 1.15 give the maxima positions, 
at f/fc = 1.05,2.05 and 3.15 for £ =2, 3 and 4 respec- 
tively. This numerical factor depends on the geometry of 
the excitation mode (bouncing in this case) : e.g. another 
(smaller) factor is obtained when the droplet is stuck on 
a vibrated solid surface [HI, [lH . 

As shown in [8^, the bouncing ability of droplets is 
due to the cooperation of (i) the droplet deformation, 
that stores potential energy, and (ii) the vertical force 
resulting from the squeezing of the intervening air layer 
between the droplet and the bath. The forced motion of 
the bath provides some energy to the droplet, a part of 
which helps the droplet to bounce (translational energy) 
while the other part increases internal motions inside the 
droplet, which are eventually dissipated by viscosity. The 
proportion of energy supplied to the translational / inter- 
nal motion varies with the forcing frequency (i.e. when 
the energy is provided in the oscillation cycle). Minima 
(resp. maxima) in the threshold curve correspond to a 
maximum (resp. minimum) of the translational to inter- 
nal energy ratio. Maxima may be related to the cut-off 
frequency recently observed and theoretically explained 
by Gilet et al. [8] . This fact is confirmed experimentally, 
since the maxima in threshold correspond to the bound- 
aries between modes. At those points, corresponding to 
a forcing frequency of 1.15 fn, the most energy is spent 
in internal motions : the droplet resonates and absorbs 
energy. 

The droplet behaviour has been observed in the vicin- 
ity of the first maximum in the threshold curve. At 115 
Hz, the droplet moves along a linear trajectory analo- 
gously as the walkers observed in [2]. In this latter case, 
50 cSt droplets were produced on a 50 cSt silicon oil 
bath. The mechanism for the walking motion is the in- 
teraction between the droplet and the bath surface wave 




FIG. 4: Rolling motion of the droplet (top) Mode Y2 of 
a bouncing droplet (z^ = 1.5 cSt, R — 0.765 mm) observed at 
/ = 115 Hz and F = 4.5 > Vth- Frames are separated by 1 ms. 
The droplet is rolling towards the left, (bottom) The rolling 
motion is revealed by small bright particles spread inside the 
droplet. 



generated by the bouncing. Such a mechanism does not 
hold in the case of a 1000 cSt silicon oil bath. Indeed, 
the generated waves are rapidly damped and cannot be 
responsible for the motion of the droplet. A movie of 
the moving droplet has been recorded using a high speed 
camera (Fig. 5). The images are separated by 1 ms. 
The deformation is not axi-symmetrical as in the reso- 
nant minima of the threshold acceleration curve. On the 
other hand, the mode is related to ^ = 2 and m = 1. 
This mode is characterized by two lines of nodes that are 
orthogonal. It results the existence of two fixed points 
located on the equator of the droplet. As the the line 
of nodes does not follow the axi-symmetrical geometry, 
those fixed points move. More precisely, they turn giv- 
ing the droplet a straight direction motion. As far as 
Rayleigh frequencies are concerned, we observe a break 
of degeneracy for the m parameter! 

Some tracers have been placed in the droplet. The re- 
fiection of the light on the particles allows to follow the 
inner fiuid motion, which clearly reveals that the droplet 
rolls over the bath surface. The initial speed v of the 
droplet has been measured in respect to the forcing pa- 
rameters (amplitude and frequency). A scaling is found 
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FIG. 5: Initial speed of the roller for various frequen- 
cies (indicated in the legend) as a function of the char- 
acteristic speed V* . 

when considering that the phenomenon can occur only 
above a cut-off frequency /o ~ 103 Hz and the amplitude 
threshold AthU) = rth{f)g/{4Try^) as 

v = 2na{A-Ath{mf-fo) (2) 

where a is a constant. The initial speed is represented 
versus the characteristic speed v"^ = {A — Ath{f)){f — /q) 
in FiglH The proportionality between the initial speed 
of the roller and is remarkable. Three characteristic 
frequencies are involved in the rolling process : (i) the 
excitation time 100 Hz, (ii) the rotation of the liquid in- 



side the droplet about 10 Hz and (iii) the translational 
motion about 0.2 Hz (the droplet travels a distance cor- 
responding to its circumference once per 5 seconds). The 
three processes are interrelated, i.e. the deformation in- 
duces the internal motion which induces the rolling and 
the translational motion. 

The deformation of low viscosity bouncing droplets is 
emphasized on a high viscosity bath : the droplet oscil- 
lations are much less damped than the bath oscillations. 
Depending on the forcing frequency, droplets need a dif- 
ferent amount of supplied energy to achieve a sustained 
periodic bouncing. When the forcing frequency corre- 
sponds to a multiple of the eigenfrequency of the bounc- 
ing droplet if ^ {i — l)/c), the supplied energy is mainly 
lost into internal motions, so the threshold acceleration 
Tth is maximum. On the other hand, Tth is minimum 
when ///, = i - 3/2. 

A new self-propelled mode has been discovered for forc- 
ing frequencies between the i = 2 and i = 3 modes 
(between roughly 100 and 150 Hz) for a sufficiently 
low droplet viscosity. This excited mode is a non-axi- 
symmetrical mode Y2 characterized by an internal ro- 
tation of the fluid inside the droplet. The droplet rolls 
over the vibrated bath ! This droplet displacement tech- 
nique is more adapted to large low- viscosity droplets; it is 
therefore of considerable interest for microfluidic applica- 
tions: manipulating aqueous mixtures without touching 
them. 
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